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I. INTRODUCTION 

In a variety of recently discovered materials, supercon- 
ductivity apparently arises directly from the electron cor- 
relations themselves. However, these materials are com- 
plex, and many material specific details can play a role 
in the mechanism of superconductivity. The problem is 
greatly simplified in the weak coupling limit, where we 
recently showed^ that an asymptotically exact treatment 
of the problem is possible, valid in cases in which the su- 
perconducting state emerges at low temperatures from a 
well-formed Fermi liquid. Nonetheless, even under these 
circumstances, the character of the superconducting state 
and the transition temperature depend in a complicated 
way on details of the band-structure, both near and far 
from the Fermi surface. 

To the extent that there are basic principles at work 
underlying the mechanism of unconventional supercon- 
ductivity, it would be a great advance to find simple 
model systems which exhibit such behavior. Here, we 
consider the possibility of unconventional superconduc- 
tivity in some model systems with particularly simple 
electronic structures, where controlled theory is possible, 
and where, conceivably, experimental tests of the theory 
are feasible. Specifically, we consider circumstances in 
which superconductivity may occur in a two dimensional 
electron gas (2DEG) in a high mobility heterostructure. 
Here, due to the the stiffness of the lattice and the lim- 
ited phase space for electron-phonon scattering, electron- 
phonon coupling is probably negligible, and the single- 
particle dynamics can be treated accurately within a ro- 
tationally invariant effective mass approximation. More- 
over, the strength of the correlations can, to a large ex- 
tent, be tuned by varying the electron density. 

The possibility of an electronic pairing mechanism in 
systems with rotational invariancc was first put forth in 
a seminal paper by Kohn and Luttinger—. Although U, 
the bare interactions among electrons arc repulsive, there 
are effective attractive interactions that arise at 0{U 2 ). 
Kohn and Luttinger focussed on the portion of the ef- 
fective attractions associated with the non-analyticities 
in x(o)i the particle-hole susceptibility, at momentum 
q = 2k f which reflect the sharpness of the Fermi surface 
at zero temperature. More generally, what is required for 



this mechanism to work is strong q dependence of x(l) 
for q < 2kp. Indeed, the Kohn Luttinger instability of 
a 3 dimensional rotationally invariant system results in 
the formation of a p-wave superconducting ground state 
due to the peak in xil) near q = 0. While this result is 
valid only in the weak-coupling regime where U <C Ep, it 
is widely believed that the p-wave ground state obtained 
this way is adiabatically connected to the more realistic 
(and more strongly correlated) example of Helium-3^. 

However, the Kohn-Luttingcr effect is exponentially 
weaker in a rotationally invariant 2DEG^, due to the fact 
x{q) is independent of momentum for momenta q < 2kp. 
It was later shown that at 0(U 3 ), the 2DEG does exhibit 
a pairing instability^. Still, at least in weak-coupling, 
electronically mediated superconductivity in the 2DEG 
is negligible. 

In this paper, we show that by perturbing the 2DEG, 
it is possible to significantly enhance the superconduct- 
ing transition temperature by engineering circumstances 
in which instabilities arise at 0(U 2 ) in perturbation the- 
ory. We present asymptotically exacfei weak coupling so- 
lutions of the superconducting instability in several sys- 
tems that are variants of the simplest, rotationally in- 
variant 2DEG. As a first example, we show that par- 
tially spin-polarizing the 2DEG produces a non-unitary 
p + ip superconductor. Y. Kagan and A. Chubukov have 
previously addressed this problem using an expansion in 
powers of the electron concentration^, and their result 
reduces to ours in the weak coupling limit. As a second 
example, we consider the 2DEG in a semiconductor het- 
erostructure quantum well with two populated subbands. 
We show that this system can possess both p-wave and 
d-wave ground states and present the phase diagram of 
this system. 

This paper is organized as follows. In the next section, 
we review the method developed in ReLi and discuss 
its straightforward generalization needed for the present 
context. In Section III, the effect of partially polarizing 
the 2DEG is studied. In Section IV, we consider two 
subbands in a 2DEG quantum well. Technical details of 
the various calculations are presented in the Appendix. 
In a forthcoming paper— we will consider a variety of 
slightly more complex situations pertinent to particular 
semiconductor heterostructures. 



II. PERTURB ATIVE RENORMALIZATION 
GROUP TREATMENT OF 
SUPERCONDUCTVITY 

In this section, we review the prescription of Refi and 
discuss its generalization to the present context. We in- 
tegrate out high energy modes in two steps. In the first 
step, we integrate out all modes outside a narrow range 
of energies Q about the Fermi energy. f2 is not a phys- 
ical energy in the problem, but rather a calculational 
device. It is chosen large enough so that the interactions 
can be treated perturbatively but small enough that it 
can be set to zero in all non-singular expressions with- 
out causing significant error, i.e. it is chosen to satisfy 
the inequalities pU 2 ^> /icxp{ — [1/pU]}, where p 

is the density of states at the Fermi energy and p is the 
Fermi energy. The effective interactions generated in the 
process then serve as the "bare" interactions in a second 
step, in which the remaining problem is solved using the 
perturbative renormalization group procedure of Shankar 
and Polchinski^. T c is, up to an unknown multiplica- 
tive constant, given by the energy scale, T*, at which a 
relevant interaction grows to be of order 1. It was shown 
by careful analysis of perturbative expressions up to 4th 
order in the interaction strength that the resulting ex- 
pression for T* is independent of Q. 

The analysis of Refi leads to the following prescrip- 
tion for computing the leading order asymptotic behav- 
ior of T c for weak interactions: First, compute the effec- 
tive interaction in the Cooper channel at energy scale f2, 
T^ a \k,k'), to second order in the interactions. Here, k 
and k' denote points on the Fermi surface, and T is the 
vertex for scattering a pair of electrons with momenta 
k and — k to states with momenta k' and —k', where if 
there are multiple band indices, the subband index is im- 
plicitly determined depending on whether the momenta 
are on one Fermi surface or the other, and where there 
is a different matrix depending on whether the electron 
pair forms a spin singlet (r' s ^) or a spin triplet (r^*'). 
We then construct the related dimensionlcss matrix 



W = p^v/v(k)r^(k,k')y/v/v(k'), (1) 
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where v(k) is the magnitude of the Fermi velocity on the 
Fermi surface of the corresponding subband, and p is the 
total density of states at the Fermi energy. Manifestly, g 
is a real, symmetric, hence Hermitian matrix, so it has a 
complete set of eigenstates and eigenvalues, 



E (os) ,(a,m) _ \(a,m)Aa,m) 
y k.k'^k' V k 



(2) 



Among all the possible solutions, we identify the most 
negative eigenvalue, 



Then, 



A = Min A (a ' m) , A < 



T c ~/zexp[-l/|A|]. 



(3) 



(4) 



III. PARTIALLY POLARIZED FERMI SURFACE 

As a first example, we consider a partially spin polar- 
ized 2DEG with short-ranged repulsive interactions: 

H = Hq+Hi 



ff ° = E/ (^vw^w^w 
i ki<?k 2 <?k 3 ^ Uj _ ^ t 



Hi — U , 

(2tt) 6 

where &4 = fci + k-2 — &3, 



^(fel)Vl(fe2)^(fe 3 )Vt(fe4) 



(5) 



(6) 



and h is a mean-field that renders the ground state 
spin-polarized. Such a partially-polarized system can 
occur in a narrow- well semiconductor heterostructure in 
the presence of a parallel magnetic field (in which case 
h = g/isH||), or in a ferromagnetic phase with spon- 
taneously broken symmetry, such as probably occurs in 
the Hubbard model away from half-filling in the strong- 
coupling limit U 3> fl 10 ' 11 . (However, in the latter case, it 
requires something of an intuitive leap to treat the resid- 
ual interactions beyond those that produce the mean- 
field h as "weak.") 

Since the Fermi surfaces are spin-polarized, singlet 
pairing is suppressed, so the leading superconducting in- 
stability will therefore be in the spin triplet channel. We 
first consider the limit in which there is no spin-orbit 
coupling, in which case, the two particle scattering am- 
plitude is a separate function for each spin-polarization. 
As derived in the appendix, 

r t (M) = -u 2 xi(k-q) 

T i (k,q)=-U 2 Xf (k-q) 

where x<? i s the contribution of spin er electrons to the 
susceptibility. 

In the case of a rotationally invariant system with 
= k 2 /2m+ah, Vf tCr (k) — kf^jm and p a = p = m/2ir 
is independent of the spin-polarization. Therefore the 
matrix ~ defined in the previous section is 



pT C7 {k,k') 



(7) 



The particle-hole susceptibility for this system has the 
following well-known form (see the appendix) : 



RevV - {2k Fa y 



(8) 



Thus, Xer(q) is a constant for q < 2kp ai has a derivative 
discontinuity at q = 2kp a , and vanishes as 1/q 2 when 
q » 2k Fa . 
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The rotational invariance of the problem implies that 
the triplet eigenfunctions are of the form 



V£" l (fc) = V(W)coe (m0- k ) 



(9) 



where m is an odd integer. The eigenvalue problem for 
this system therefore reduces to the integral expressions: 

A (m ' t} = - P U 2 J Jl Xi (2fc Ft |sin (0/2)|) cos (m0) 

/rlf) 
—Xt (2fc n |sin(0/2)|) cos (m0) 

(10) 

where 9 is the angle relative between k and q. 

Without loss of generality, we suppose that k F i < fc*t- 
For any k 7 q on the smaller (spin-down) Fermi surface, k— 
q < 2k F i < 2k Ff so the effective interaction, ~ X^{k — q), 
is a constant. Therefore, it follows that A m j, = for 
all m or in other words, the smaller Fermi surface has 
no superconducting instability to 0(U 2 ). ( Presumably, 
\ m ,i~0(U 3 ). ) 

Conversely, the effective interaction between electrons 
on the larger (spin-up) Fermi surface is ~ Xl(k 0), which 
does depend on the relative position of the incoming and 
outgoing electrons on the Fermi surface. Using Eq. [SJ the 
above expression for the eigenvalue on the larger Fermi 
surface becomes 



AmfW) = / dO- 

ir Ja„ sin ; 



cos (me) (11) 



where rj — (kpi/kp-f), < rj < 1, and 
can be seen from the equation above, A TO -f(0) 



2 sin 1 rj. As 
A m f(l) 

!). That is, in the limit where the Fermi surface is cither 
completely polarized, or completely unpolarized, there 
is no superconducting instability to 0(U 2 ). However, 
for intermediate values of the polarization, the integral 
above yields 



(12) 



which is clearly negative for all intermediate values of rj. 
This is the main result of this section: by polarizing the 
Fermi surfaces in two dimensions, there is a significant 
enhancement of p-wave superconductivity. The optimal 
pairing strength occurs when rj = 1/2, so that 

Max [A m -f (?/)] = Aitfa = 0.5) = (13) 

( Note that Eq. 7 of Ref£ reduces to this result in the 
limit of weak interaction. ) For completeness, we quote 
the next leading eigenvalue, which corresponds to the f- 
wavc (i.e. m = 3) solution: 

ht(v) - ~P 2 U 2 V [1 - r, (3 - 4rf + 2r? 4 )] (14) 

which is not symmetric about the point rj = 0.5. 



Weak but non- vanishing spin orbit coupling will gener- 
ically change this situation, since superconductivity will 
be induced in the minority fluid by the proximity effect as 
soon as the majority fluid becomes superconducting. In 
2D, this induced superconductivity will generally track 
the fundamental order parameter. 



IV. TWO SUBBANDS IN A 2DEG 

In this section, we consider the case of a 2DEG in 
a semiconductor heterostructure having two subbands, 
with Hamiltonian: 



H = Ho+m 

fl » = EE 



d 2 k 

[2ftf 



#1 - ^2^2 V a£d 
a..d ct.cf' 



d 2 k 1 d 2 h 2 d 2 k 3 
J2ftf ^ 



(15) 



where a is the subband index and is used to distinguish 
the smaller (a = 1) and larger (a = 2) Fermi surface, and 
ffc.a = k 2 /2m + 6 a with S± = and 82 > 0. The inter- 
actions are assumed to be short-ranged, consisting of an 
intra-band repulsion U, an inter-band repulsion V, and 
an inter-band pair-hopping amplitude J. The interaction 
matrix in the basis (lalcr', la2a\ 2<j\o' 2a2a') is thus 



v: 



b.cd 



where 





/ u aa > 








J crcr' 




















V aa < 










V J (7<7 ! 
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J <7(7 ! 


■7(1- 


- <W 


) 



(16) 



V„„. = V 



(17) 



As before, rotational invariance enables us to label the 
eigenstates by the eigenvalue of the rotation operator, 



(18) 



where the complex amplitude (f> a depends only on the 
subband index associated with k, is the angle between 
k and an arbitrarily defined x axis, and m must be an 
even integer in the singlet channel and an odd integer 
in the triplet channel. Consequently, for each integer m, 
rotational symmetry reduces the eigenvalue problem to a 
2x2 problem, 



E 



fa.a'S 



(m) 



■) 



(19) 
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where 



where U a t are momentum independent interactions. 



~(m) _ 
9 a J = 



_ f qfa) -im» irafl 

2tt L 2?r M' 



(20) 



where y = s (singlet) for m even and y = t (triplet) for 
m odd. The most negative eigenvalue for fixed to is 



A ( ' 



' ~(m) 
Sl.l 



-(mV 
.92,2 



(21) 



'-(m) 



-(to)' 
52,2 



l~ « 



We first consider the spin triplet channel (to odd) 
which is only a slight extension of the result obtained 
for a partially polarized Fermi surface. As shown in the 
Appendix, for odd to, the effective interaction is diagonal 
in the subband index, and depends on U, V , but not J: 

<?S(M) = -pU 2 Xl , 1 (k-q)-2pV 2 X2a {k-q) 
ffS(M) = -pU 2 X2, 2 (k-q)-2pV 2 Xi,iCk-q) 



91,2$, q) 

where 



Xa,b(k) 



d 2 p f(e p +k,a) - /(e P ,b) 

(27T) 2 



(22) 



(23) 



is the particle-hole susceptibility generalized to the two 
band system. The intraband susceptibilities are precisely 
the same functions used before: 



Xa,a(q) 2 



ReV? 2 - (2k Fa y 



(24) 



with the subband index playing the role that the spin 
played in the previous section. Therefore, we may sim- 
ply transcribe the results found in the previous section to 
the present context. The band which forms the smaller 
Fermi surface (a = 1) does not exhibit a superconducting 
instability to 0(U 2 ). The larger Fermi surface exhibits a 
triplet p-wave instability with a pairing strength deter- 
mined solely by V: 



A (1) (ry) = -Ap 2 V 2 r}{l-r]) 



(25) 



where r\ = (A1.Fi/fcF2)- 

In the spin-singlet channel, the matrix g has off- 
diagonal components: 



g[ s }(k,q) = pU l ~2pV 2 X 2,2(k-q) 
9 { 2%{kq) = pU 2 ~2pV 2 X i,i{k-q) 

9i.kk,q) = pU 12 + P vj \xidk + q) + xiA k -q) 



(26) 



t/i = U + U 2 P 1 (n) + J 2 P 2 {Q) + U 2 xi,i(k + q) 

u 2 = u + u 2 p 2 (n) + j 2 p 1 (n) + u 2 x2Ak + q) 

U u = UJ[P 1 (fl)+P 2 (n)} (27) 
where the particle-particle susceptibility, 



d 2 q 2/(e,,„)-l 



plog 



(2tt) 2 tfl-2e q , a 
E F -S a 



9. 



(28) 



is a momentum-independent constant which diverges log- 
arithmically at low energies. Despite this divergence, the 
second order contributions to U a are unimportant, since 
they do not enter the gap equation for any m 7^ 0, and 
any s-wave solution is already killed by the first order 
terms proportional to U. 

For to > and even, the effective intra-band interac- 
tion depends only on the interaction V and is non-zero 
only for the larger Fermi surface, whereas the inter-band 
interaction depends both on V and J: 



-(to) 

01,1 

~(m) 
92.2 



-(to) 
9l,2 



ko 







-2V 2 p J I^Xm (2fcF 2 |sin (0/2)|) cos (m0) 



d6 

2V 'J 'P I — xi,2 (h) cos (m9) 

Z7T 



<F2 



y / (l-7 ? ) 2 +47 ? sin 2 (0/2). 



(29) 



The explicit expression for the inter-band susceptibility 
Xi,2{q) is derived in the Appendix. Since the to = 
eigenvalues are always positive, the dominant singlet in- 
stability is in the d-wave (to = 2) channel. The quantity 
g 2 2 is obtained by computing 



, n s r dO Re\/sin 2 % - rf- 

gS 2 \ = -2V 2 p 2 \ V , I cos (20) 

J- 2w sin g 



-V 2 p 2 rj{'q-l){rf +r)-l) 



(30) 



The intcrband interaction g\ 2 is also obtained using Eq. 



-(2) 
01,2 



where, for < x < 1, 



VJp 2 
2tt 



$(77) 



2 sin 1 x 



(31) 



(32) 



This function is discussed in detail in the Appendix. 
An important property of $(x) is that it is a monotoni- 
cally increasing function of x for < x < 1 (see Fig. [3]). 
Therefore, the effective intcrband scattering grows with 
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r/. The pairing strength in the d-wavc channel is obtained 
from these quantities via 

A< a >fo) = ^-iVW+*W (33) 

Having derived closed form expressions for the p-wave 
and d-wave pairing strengths, we can construct the phase 
diagram, shown in Fig. [TJ The phases are labeled ac- 
cording to the symmetry of the most negative eigen- 
value, so the phase boundaries are the lines at which 

A (2) = A (l) < 

Since the d-wave and p-wave eigenval- 
ues are both negative for all < 77 < 1 , where one phase is 
stable, the other is metastable. It would require different 
methods of analysis to completely characterize the phase 
competition. However, in weak coupling, the phase with 
the larger |A m | has an exponentially larger T c , and so 
gaps the entire Fermi surface at temperatures far above 
the putative transition temperature of the subdominant 
phase. Thus, a BCS mean-field treatment of this prob- 
lem would suggest that at low temperatures, there is a 
direct, first order transition, or at most an exponentially 
narrow region of phase coexistence between the extremal 
pure d-wave and p-wave phases. 

Since the pair-hopping term only affects spin-singlet 

superconductivity, and since > \gi 2 l\, it follows that 
for J = 0, the p-wave solution always remains the favored 
ground state, as can be seen from Fig. [TJa). However, as 
the interband scattering is enhanced, the d-wave pairing 
strength grows. Since the interband scattering increases 
monotonically as a function of both 77 and J, it is seen 
that for sufficiently large values of either parameter the 
p-wave superconductivity gives way to a d-wave ground 
state. In Figs. HJc-d), we show how the magnitude of 
A 1 - 1 ) and A^ 2 ) depend on 77, from which one can see that 
T c is maximal in the p-wave channel when ?/ = 0.5. How- 
ever, when J =/= 0, the d-wave channel grows monotoni- 
cally with r\ and ultimately overtakes the p-wave pairing 
strength as 77 — > 1. Note, however, that 77 can never 
equal unity in this context, since it is determined by the 
thickness of the quantum well. 

V. DISCUSSION 

We have obtained analytical expressions for various 
unconventional superconducting ground states of a clean 
2DEG in the presence of weak, short-ranged repulsive in- 
teractions. Ultimately, to make contact with experiments 
involving real 2DEGs, we must take into account the 
Coulomb interactions. In the small r s limit, the Coulomb 
interactions are sufficiently well screened that it may be 
reasonable to treat them as weak and short-ranged^. 
We thus imagine we can relate the physical problem to 
a problem with short-ranged interactions and speculate 
on two ways in which unconventional superconductivity 
could be found in the 2DEG in physically realizable semi- 
conductor heterostructurcs. (We shall present more com- 
plicated examples in a forthcoming publication^) 



In the first scenario, an in-plane magnetic field is ap- 
plied to partially polarize the 2DEG in a narrow quantum 
well. This system is predicted to exhibit p-wave pairing 
with a transition temperature which is non-monotonic in 
the magnetic field. The optimal transition temperature 
is obtained for a magnetic field at which the ratio of the 
distinct spin Fermi momenta is 77 = 1/2. In the sec- 
ond scenario, the 2DEG is confined to a relatively broad 
quantum well, and the density is tuned to the range in 
which two transverse subbands are occupied. For fixed 
total electron density, the ratio, r/ 2 , of densities in the 
two subbands increases with increasing thickness w of 
the quantum well. When this ratio is small, a p-wave 
groundstate arises, with a T c that rises sharply with in- 
creasing 77 so long as 77 < 1/2. However, this gives rise to 
a d-wave ground states above a certain critical thickness. 

Insight into the dependence of V/J on the thickness, w, 
is obtained by considering the Coulomb interactions. A 
simple estimate shows that for kpw <C 1, V ~ e 2 /kp and 
J ~ V (wkp). Therefore, for thicker quantum wells, J 
becomes increasingly important and favors d-wave pair- 
ing whereas thinner quantum wells should exhibit p-wave 
pairing. Depending on the ratio of V/J, the optimal T c 
occurs either for 77 m 1/2 (p-wave) or for the largest pos- 
sible 77 (d-wave). In both cases, T c ~ Ep exp [—a/ (pV) 2 ] 
where a is an 0(1) constant. We have found that for d- 
wave superconductivity in the 2 subband system, values 
as low as a ~ 1 are within reach. 

Three practical considerations warrant mention. Due 
to the unconventional nature of the superconductvity, it 
is very fragile to even weak quenched disorder. There- 
fore, the results presented here arc likely to be realized 
only in the purest samples with mean free paths exceed- 
ing the Fermi wavelength by several orders of magnitude. 
Furthermore, for small r s , the plasma frequency is small 
compared to the Fermi energy, i.e. u p ~ y/r^Ep, so 
even if it is reasonable to treat the interactions as short- 
ranged at low energies, this approximation is certainly 
not valid all the way to the Fermi energy. Finally, since 
the transition temperatures are exponentially low in the 
effective interactions, ultimately the superconductivity 
studied here is likely to be observable only in the regime 
r s ~ 1, where the long-range character of the Coulomb 
interaction may not be negligible, and where, even for 
short-range interactions, a well-controlled solution to the 
problem is unfeasible. We therefore are forced rely on the 
hope that the asymptotic results smoothly extrapolate to 
the intermediate coupling regime, where it is conceivable 
that these states can be observed in experiment. 

With these caveats, we turn to the most uncertain part 
of the discussion, and make the following crude quan- 
titative estimate of T c based on our calculations: We 
identify V with the Fourier transform of the Coulomb 
interaction evaluated at kp, i.e. V w e 2 ir/kp, from 
which it follows that pV w (r s /4). Since we are go- 
ing to extrapolate to r s ~ 1 in any case, we simply 
ignore subtleties associated with the small value of oj p . 
Then, T c ~ Ep cxp[— a(4/r s ) 2 ], where, for optimal cir- 



ri 




FIG. 1: Phase diagram of a 2DEG having two subbands. (a) Phase diagram for fixed 77 = fc.pi/fcF2 = 0.5 as a function of 
the dimensionless couplings pV and pJ. U does not enter the problem except in that it is responsible for the suppression of 
s-wave pairing, (b) Phase diagram for fixed pV = 0.1 as a function of rj and pJ. c) The dimensionless strength of the pairing 
interaction in the p-wave (solid line) and d-wave (dashed line) channels for fixed J = 0. d) Same as c), but for J = 0.2V. 



cumstances awl. 
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Appendix: Perturbation Theory 



In this section, we derive the perturbative expansion 
of the effective interaction in the Cooper channel for the 
problems studied in the main text. Due to the presence 
of multiple bands, adopting a more compact notation 
enables us to treat all of the above problems in a unified 
fashion. Wc consider a Hamiltonian of the form 



(2c) 



(2d) 



(2e) 



H = H +H 1 

H = ^2ek,acl a C k ,a 
k.a 

H l = Yl Yl V ab,cdcl i a cl 2 b C k3tC C kiA (A.l) 
h\ ,&2 ,^3 a,b,c,d 

where /C4 = k\ + &2 — k%. The Latin subscripts denote 
a collective set of "band indices" which label the energy 
cigenstatcs. In the problem of partially polarized Fermi 
surfaces, they simply label the majority and minority 
spin bands. In the problem of multiple subbands in a 
quantum well, it indexes the subbands, and in the prob- 
lem involving Rashba spin-orbit coupling, the Latin index 
refers to the positive and negative helicity subband. The 
bare interaction Hi is to be interpreted as a matrix; it's 



FIG. 2: Diagrams which contribute to V(k,q), shown to 
quadratic order in the interactions. Each of the incoming 
quasiparticles has momentum ±fc and band index a. The 
outgoing electrons have momentum ±g, and are in band b 
(with the exception of the first diagram, the momenta and 
band indices of each diagram are not shown). Both intra- 
band and inter-band scattering processes contribute to the 
effective interaction. 



rows labels the outgoing states and its columns label the 
incoming states. 

Higher order scattering processes are derived using di- 
agrammatic perturbation theory in the usual manner—. 
In addition to integrating over the internal momenta, the 
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band indices of any internal line are also summed over, 
weighted by the appropriate component of the interaction 
vertex, as will be made clear from the examples below. 

The primary quantity of interest here is the two- 
particle scattering amplitude in the Cooper channel, de- 
noted r(fc, q), which is the amplitude for scattering a pair 
of electrons with momenta ±fc into a pair with momenta 
±q. If the system at hand possesses inversion symme- 
try (so that the Kinetic energy consists of terms that are 
even powers of momentum), superconducting states can 
be classified has being even or odd parity states; the for- 
mer include (e.g. s-wave, d-wave, etc.), and the latter 
include, (p-wave, f-wave, etc.) instabilities are perfectly 
decoupled from one another. If, in addition to inversion 
symmetry, spin-rotation symmetry is also preserved, then 
the scattering amplitudes in the singlet channel consist of 
processes in which the incoming electrons have opposite 
spin polarizations, whereas in the triplet channel, they 
have identical spin polarizations. On the other hand, 
when inversion symmetry is broken by the presence of 
Rashba spin-orbit coupling, there is no sharp distinction 
between even and odd parity pairing. However, since 
the Rashba coupling breaks the 2-fold degeneracy of sin- 
gle particle states at each momenta, there is pairing be- 
tween states of opposite helicity (i.e. opposite momenta 
and opposite in-plane component of the spin). 

Figure [5] shows the lowest order Feynman diagrams 
which contribute to T(k,q). Generally, all of these dia- 
grams contribute both in the singlet and triplet channel. 
The diagrams are each equivalent to 

1 : Vbb,aa 

2a: ^ Vbb,ccVcc, aa / G c (-p)G c (p) 

c Jv 

2b : V bc,adVdb,ca / G c ( P )G d { P + k + q) 

c,d J P 

2c: ^ Vbc,ad V db,ac / G c (p)G d (p + q - k) 
c.d J P 

2d : Vbd,caV bc ,ad / G c (p)G d (p + k-q) 

c,d J p 

2e: - ^ V bd ,caV c bMd / G c (p)G d (p + k - q) 

_ J J V 



c . d 



where 



and 



duj p d 2 p 
(2^) 3 



G( P ) 



(A.2) 
(A.3) 

(A.4) 



is the single particle Green function of the non- 
interacting system. We next apply this general formalism 
to the three problems studied in this paper. 



1. Partially polarized Fermi surfaces 

Although this problem is rather simple, and the Feyn- 
man rules for a single band system are sufficient, we will 
apply the notation above to this problem. This will cer- 
tainly prove to be valuable in the more non-trivial ex- 
amples studied thereafter. The interaction vertex, in the 
basis (ttjti)4T)44-)) is 



b.cd 



/ 
U 
U 

Vo o o o 



(A.5) 



Note that in this problem, there are no processes which 
scatter two electrons from one Fermi surface, to two elec- 
trons in a different Fermi surface. Therefore, T(k,q) 
is diagonal in the band index - which is just a long- 
winded way of saying that only equal spin pairing (i.e. 
spin-triplet pairing) can occur. Indeed, the only diagram 
which contributes to the effective interaction is 2c, which 
yields: 



T n (k,q) 

r u (M) 



u 2 Xl (k 
-u 2 Xt (k 



q) 

q) 



o(n ) 
o(n ) 



(A.6) 



where 



X „(k) = I G a {p)G a {p + k) 

d 2 p /(ep+fc.q) - f{e p , a ) , , 

(ZTT) €p+k,a - (p,a 

is the non-interacting susceptibility of each spin band. 

2. Multiple subbands in a 2DEG quantum well 

For the problem involving two-subbands in 
a quantum well, we choose the basis to be 

(lcrl(T',lcr2cr',2crlcr',2cr2cr') and 



V, 



/ u aa > 








J (7(7 








Voo> 








v aa , 








V Jaa' 











(A.8) 



where 



J 70 1 

V aa > 



J(i-<W) 
v 



(A.9) 



In this basis a and a' refer to the spins of the incom- 
ing electron states. Since there is inversion and spin- 
rotational symmetry in this problem, we can study the 
effective interaction in the singlet (a = —cr') and triplet 
(cr = a') channels separately. We shall refer to the ef- 
fective interaction as r s m where the subscript stands for 
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singlet and triplet. Having specified the spin polariza- 
tions of the electrons, r s ( t ) will still be a 2 x 2 matrix 
due to the presence of two subbands. We will denote this 
as 

r (h n \ - ( r s(t)(ki,qi) r s(t) (fei,q 2 ) \ /a i rt\ 

with the understanding that fc a denote momentum states 
associated with band a. 

Next, we state the contributions from each of the dia- 
grams in Fig. [2J First, in the singlet channel, 



r.(i) = 

r s (2a) = 

r s (26) = 

r s (2c) = 

r s (2d) = 

r s (2e) = 



U J 
J U 

U 2 P 1 + J 2 P 2 UJ(Pi+P 2 ) 
UJ{P 1 +P 2 ) U 2 P 2 + J 2 P 1 

U 2 X i,i{ki+qi) VJ X i. 2 {ki+q 2 ) 
VJx2,i{k 2 + qi) U 2 x 2 , 2 ( k 2 + 92) 

VJ X i. 2 {ki -q 2 ) 

qi) 

vJxiA k i ~ 92) 

qi) 



^X2,i(fc 2 


^X2,l(fe 2 

2V 2 X 2, 2 (ki - qi 






2V 2 X is(k2-q 2 ) 



(A.ll) 



where 



Pn 



Xa,b( k ) 



G a ( P )G a (-p) 

Pa log [A/n Q ] + o(n Q ) 

G a ( P + k)G b ( P ) 
d 2 P /(e P +fe,a) ~ /(e P ,h) 

(27r) 2 e p +fc, a — £p,fc 



(A.12) 



are the particle-particle and particle-hole susceptibilities, 
respectively. 

In the triplet channel, only diagram 2e has a contribu- 
tion and I\(fe a , qb) is diagonal in subband index: 

r t (fei,qi) = -U 2 xi,i{ki-qi)-2V 2 X2, 2 {ki-qi) 



r*(fc 2 ,g 2 ) 
r t (fe 1 ,q 2 ) = 



-£^ 2 X2,2(fc2 - q 2 ) - 2V 2 xiA k 2 - 92) 

(A.13) 



Having computed the effective interaction T s , t , we define 
the quantity 



g s ,t{k a ,qb) 



Vf(k a ) 



r(fc Q ,g b ) 



Vf(qb) 



(A.14) 



which is also a matrix whose row and column indices are 
the set of momentum states on the Fermi surface. 
3. Some relevant integrals 

In this section, we compute the particle-hole suscepti- 
bility matrix of the two subband problem: 



Xa,b{k) 

We let 



d 2 P f((- P +k,a) - f(e P ,b) 



(27T) 



(A.I5) 



£fc,2 = (k + A 



(A.16) 



and set A = {k F1 — kp 2 )/2m > without loss of gener- 
ality. It follows that the intraband susceptibilities are 



Xaa{q) = 2 



2jt 



dO 



kdk 



(2n) 2 J e q + kqcos8/m 



(A.17) 



where Ufi = (2mfi) 1 ^ 2 and k F2 = (2m(/x + A)) 1 ^ 2 . The 
integrals are standard, resulting in the following: 



m 
2^ 



1 - 



Rcv/q 2 - {2k Fa ) 



(A.I8) 



The interband susceptibility is 

Xl,2(<?) = 



kF1 kdk r2 * 



o (27r) 2 J e q — A + kqcos9/m 

kF2 kdk r 271 do 



(27r) 2 Jq e q + A + kq cos 9 /m 

(A.I9) 



These integrals are straightforward and they produce the 
final result: 



Xl,2(<?) 



ReJq 2 [I - X(q)} 2 - (2 V k F2 ) 2 ReJq 2 [1 + X(q)} 2 - {2k F2 ) 2 



[1 + A( 9 )]-|1-A(«)| 



(A.20) 
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where 



V = 



&F2 



m = *f (i - v 2 ) 



(A.21) 



The effective inter-band attraction g^' 2 in the singlet 
channel is related to this susceptibility, as discussed in 
section IIVI 



~(s,m) OT/ T 

g\ 2 = 2VJp 



dO 
2t7 



X1.2 (h) cos (m6) 



ko 



k F2 \/(l-vf +4ijsm 2 (9/2) (A.22) 



It is easy to show that the first two terms in Eq. IA.20I 
do not contribute to Xi,2{kg) since 



Jkj [1 - \{k e )f - (2r,k F2 f 



2t] I — 

= — V-sin^ 

ka 



6 (A.23) 



where Iq = kg/kp 2 . The integration in Eq. IA.22I is per- 
formed in the complex plane defining z = —id. We find, 
for the d-wave case (m = 2), 



9i, 2 = —^(V) 



2tt 



(A.26) 



where, for < x < 1, 



2 sin 



(A.27) 



and $(1) = 0. The function $(x) is shown in Fig. [3] 
It should be noted that is discontinuous at x = 1, 

which results from the singluar behavior of the interband 
susceptibility in the limit where the bands become degen- 
erate: 



is purely imaginary for < r\ < 1. Thus, 

XiM = § [1 + A (M + I 1 - W|] (A. 24) 

This in turn can be rewritten as 

C 1, cos# < r\ 
XrM= P \ M cos0>?? (A.25) 



lim limxi,2(g) = p— 

77— s-l g— >0 1 



(A.28) 



However, this feature does not have a physical conse- 
quence since there is always a non-zero splitting between 
the bands caused by the finite thickness of the semicon- 
ductor heterostructure. 
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FIG. 3: The scaling function &(rj) which determines the ef- 
fective interband interaction g s {\. 



